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ABSTRACT 


Cole  and  Huth  [1]  have  studied  the  effect  of  a  line 
load  moving  with  constant  velocity  V  along  the  surface  of 
an  elastic  half-space.  The  present  paper  treats  the  equiva¬ 
lent  problem  for  a  viscoelastic  (standard  solid)  material 
when  the  velocity  V  is  less  than  the  velocity  c^  of 
shear  waves  of  high  frequencies  {q  ->  oo)  . 

In  Reference  [1] no  solution  could  be  obtained  when  the 
velocity  V  was  equal  to  the  velocity  of  Rayleigh  waves. 

The  present  analysis  yields  a  solution  for  this  special  ve¬ 
locity.  It  also  permits  an  evaluation  of  the  effects  of 
viscosity  at  other  values  of  V  .  In  certain  ranges  these 
effects  are  minor j  but  in  other  ranges  major  differences  oc¬ 
cur  due  to  focusing  phenomena. 
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THE  EFFECT  OF  A  MOVING  LOAD  ON  A  VISCOELASTIC  HALF-SPACE 


I ,  Introduction 

The  response  of  an  elastic  half-space  due  to  loads  moving 
on  the  surface  has  been  considered  by  Sneddon  [1^  2.,  3]  ’'^ti.o 
obtained  formal  solutions  for  general  loads j  and  closed  solu¬ 
tions  for  certain  specific  cases.  Subsequently,  Cole  and 
Huth  [4]  obtained  by  other  methods  closed  form  solutions  for 
the  plane  problem  of  a  line  load  progressing  with  constant 
velocity  on  the  surface  of  the  half-space.  Miles  [5]  has  con¬ 
sidered  the  case  of  loads  with  axially  symmetric  distributions 
acting  over  a  circular  area,  the  radius  of  which  expands  with 
velocity  V,  Figure  1.  He  has  demonstrated  that  the  two-dimen¬ 
sional  Cole  and  Huth  problem,  with  an  exception,  is  identical 
with  an  asymptotic  solution  of  the  three  dijmensional  problem 
valid  in  the  vicinity  of  the  load  front  [5].  This  confirms 
the  intuitive  conclusion  that  plane  steady-state  solutions, 
which  are  relatively  easy  to  obtain,  can  be  used  (with  limi¬ 
tations)  as  approximations  in  three  dimensional  situations 
such  as  that  shown  in  Figure  1. 

The  results  of  [4]  indicate  that  the  character  of  the 

solution  depends  on  the  relative  values  of  the  velocity  V 

and  of  the  velocities  of  dilatational  and  shear  waves  in  the 

medium,  c  and  c  ,  respectively.  The  three  resulting 
p  s 

cases  were  designated  in  [4]  as  subsonic,  trans-sonic  and 
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supersonic.  Subsequently,  it  was  noted  in  [5]  that  the  steady- 
state  solution  of  [4]  breaks  down  if  V  =  the  velocity 

Xv 

of  Rayleigh  waves,  such  that  in  this  exceptional  case  the  three 
dimensional  problem  cannot  be  approximated  by  the  plane  solu¬ 
tion. 

Lorsch  and  Freudenthal  [6]  have  considered  a  related 
problem,  the  quasi-static  problem  of  a  line  load  moving  with 
constant  velocity  on  the  surface  of  a  viscoelastic  (Maxwell 
solid)  half-space. 

The  plane  steady-state  problem  of  a  line  load  progressing 

with  constant  velocity  V  on  the  surface  of  a  viscoelastic 

(standard  solid)  half-space  has  been  treated  by  Sackman  [7]j 

but  only  for  the  supersonic  case,  i.e,,  when  V  is  larger 

than  all  wave  velocities  in  the  medium.  In  the  present  paper, 

an  alternative  case  will  be  considered  where  V  is  smaller 

than  the  limiting  velocity  c  of  waves  of  high  (lim  -^  <») 

s 

frequencies . 

In  a  standard  solid  the  wave  velocities  decrease  with 
decreasing  frequency  Q,  the  limiting  velocities  for  0 

being  ^  ^  present  case  will  therefore  in¬ 

clude  subcases  depending  on  the  value  of  V  relative  to  c* 
and  c*.  The  solutions  obtained  permit  an  insight  into  the 
deviations  between  the  response  in  elastic  and  viscoelastic 
situations  in  general;  further,  a  solution  is  presented  for 
the  special  case,  V  =  c„,  where  no  steady  state  solution 
for  the  elastic  half-space  exists. 
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A  formal  solution  of  the  problem  will  be  obtained  by  in¬ 
tegral  transform  methods.  Due  to  the  complexities  introduced 
by  the  viscoelastic  behavior,  closed  solutions,  do,  in  general, 
not  exist,  and  recourse  will  be  taken  to  asymptotic  methods 
of  evaluation.  The  quantities  of  interest,  the  stresses  and 
accelerations,  will  be  obtained.  Additional  quantities,  such 
as  velocities  or  displacements,  could  be  obtained  if  desired, 
except  for  a  free  and  undeterminable  constant  in  the  displace¬ 
ments  .  The  indeterminacy  of  the  displacements  is  not  unex¬ 
pected;  it  already  exists  in  the  two-dimensional  elastic  prob¬ 
lem  [4],  and  even  in  the  equivalent  static  one  [8].  (The  mat¬ 
ter  is  discussed  in  [5].) 
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II.  Formulation  of  the  Problem 


Let  Xj  Y,  z  designate  a  stationary  coordinate  system 
in  the  half-space  z  ^  0,  Fig.  2.  A  line  load  P  of  unit 
intensity^  positive  if  dovmward^  moving  in  the  negative  x  - 
direction  with  uniform  velocity  V,  may  be  described  by 


P(x,  t)  =  6(x  +  Vt)  (1) 

where  the  symbol  6  indicates  Dirac's  function.  The  equa¬ 
tions  of  motion  are 


a 


ij,  j 


(2) 


The  stress-displacement  relation  for  a  homogeneous  iso¬ 
tropic  material,  elastic  in  bulk,  and  viscoelastic  in  shear  is 


=  (k  - 


a(u .  .  +  u .  . ) 

i,D 


(3) 


where  k  is  the  elastic  bulk  modulus.  For  the  standard  solid 
]1  is  the  operator 


M- 


(J- 


m  +  T 


1  +  T 


(4) 


In  this  relation  p,  is  the  "unrelaxed"  and  m(j,  the  "relaxed" 
shear  modulus,  while  T  is  the  relaxation  time,  (See  Ref.  9) • 
The  displacements  may  be  written  according  to  the  Helm¬ 
holtz  resolution. 


u . 
1 


+  e 


ijt 


(5) 


as  functions  of  the  potentials  $  and  In  the  present 
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and  the  provision  that 

^0  as  z  00 

The  problem  to  be  solved  is  now  specified  by  Eqs.(9)i 
(10)  and  (11). 


(11) 


III.  Formal  Solution 


The  solution  of  Eqs.  (9)>  (10)  and  (ll)  may  be  obtained 

formally  by  the  use  of  an  integral  transform  technique.  The 

steady  state  disturbance  from  loads  moving  with  a  velocity 

V  <(  c  is  expected  to  r'xtend  to  infinity  in  both  directions 
s 

such  that  the  Fourier  transform  is  appropriate. 

The  Fourier  integral  relations 

_  r  io)^ 

F(cd)  =  F(^)e  de  (12) 


—  CO 


apply  provided  F(x)  satisfies  certain  conditions  [10,  11]. 

As  a  first  step  the  Fourier  transforms  for  the  stresses 
and  accelerations  will  be  obtained. 

X  z 

Introducing  dimensionless  coordinates  ^  ,  T]  =  —  , 

Eqs.  (9)  give: 

fl  -  ]  $  =  0 

r)T|  I  L  n-ico  / 

(14) 

-  03^  (1  -  m£  ?  =  0 

■’  irr]  \  T?  m-icD  / 


while  Eqs.  (lO)  give: 

-“"K  -  2  feiS)  0)  -  2ic.  0). 

-2i(D  0)  -  (2  -  0)  =  0 


VT 

U- 


(15) 
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,2  .  Vf 


where  the  following  abbreviations  are  used: 

P  P 

"  K  *^3  ^  “^s  "  n  =  1  -  4/3K(1  -  m)  . 


The  symbols  c  and  c  indicate  the  velocities  of  irrota- 
p  s 

tional  and  equivoliominal  waves,  respectively,  in  the  limit  of 
infinite  frequency,  ^  In  addition,  Eq.  (ll)  requires 


that 


<5,^  are  to  remain  finite  as  Tj  ->  oo  . 
The  solutions  of  Eqs .  (l4)  are 


v^caA]  1-  M: 


2  l-i<X) 
L  n-ico 


.L  w?  l-io) 

"l 


+  A  e 
2 


+  B  e 
2 


To  make  the  quantities  CD  A  1-  and  cd 

u  L  n-icD 

single  valued,  branch  cuts  and  argument  restrictions  are  made 
below.  Further  analysis  requires  separate  consideration  of 
three  cases,  depending  on  the  value  of  the  velocity  V.  De¬ 
noting  the  velocities  of  irrotational  and  equivoluminal  waves 
in  the  limit  0  by  c*  =  sfnc  and  c*  = '/me  ,  the  fol- 

lowing  three  possibilities  must  be  considered  separately. 


Case 

a : 

V  < 

c* 

s 

Case 

b: 

=s< 

v< 

c* 

p 

Case 

c : 

C*  \ 

P 

v< 

c 

s 
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Case  a 


If  the  branch  cuts  shown  in  Fig.  3  used  with  the  ar¬ 

gument  restrictions 


- "  <  “9(1-  “?  ra)  <  " 

(18) 

-  7r<  arg(l-  m£  ^ )  <  tt  (b) 


then 


± 

1  (2  -  Mj)e 

p-CU^AQ 

(19) 

f  =  1  2ip,  e 

licD^AQ  ^ 


where  the  upper  signs  are  to  be  used  for  R(co)  <(  0^  the  lower 
signs  for  R((jo)  0^  and  the  following  abbreviations  are 
employed : 


Pl  =  1- 

=  1- 

l-ico 

L 

L 

n-ioo 

P?  =  1- 

—2 

=  1- 

1-icu 
m-  ioo 

A  =  (2  - 


(20) 


The  exponents  thus  chosen  have  negative  real  parts  in  the  range 
of  CO  where  they  are  to  be  used. 

Substitution  of  Eq.  (5)  into  Eq.  (3)  makes  it  evident  that 
the  stresses  may  be  written 
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ia>(2  -  m5)  i  ''“Pl 

- zrn —  ® 

IjlT^A  q 


M-T^A  q 


± 

e 


a3p^(2  -  M^) 
IJ-T^A  Q 


e 


2mp,  i  '‘“^T 

Xi 

_ _ ^ 

IjlT^A  Q 


(23) 


Case  b 


If  the  branch  cuts  shown  in  Fig.  4  are  used  with  the 
argument  restrictions  (l8b)  and 


0  <  arg(l- 


l-ico 
m-  io) 


)  <  27T 


(24) 


then 


<!> 


VT 

M-Cjd^Q 


(2  -  M^)e _ 

(2  -  mJ)2  +  40^0^ 


f  =  + 


VT 

1X00^  Q 


2ip^e 


(2  -  +  4p^p^ 


(25) 


Again,  the  upper  signs  apply  for  R((X))  “C  0,  the  lower  signs 
for  R(<x))  )>  0  and  the  exponents  have  negative  real  parts. 
Using  the  definition  (21)  and 

A-  =  (2  -  mJ)^  +  4PlP^  (26) 
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The  transforms  of  the  accelerations  are 


CM  |4 
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IV.  Asymptotic  Approximations  for  Large  Values  of  t] 

According  to  Eqs.  (13)^  (19)^  (25)  {29),  the  stresses 

and  accelerations  are,  in  general,  given  formally  by  expres¬ 
sions  of  the  form 

o  T^oah_(cjD)  00  'qcoh_^(u:.) 

I  =  J  dcjD  +  J  9^(oi)e  dto  .  (32) 

—  00  o 


Since  the  exponents  have  been  chosen  with  negative  real  parts 

and  the  functions  g(a3)  are  finite  along  the  real  axis  with 

the  exception  of  the  point  |  od  j  =  oo  when  V  =  c  ,  and  of 

R 

the  point  od  =  0  when  V  assumes  the  special  value  c*  dis- 
cussed  following  Eqs.  (39) >  the  integrals  of  Eq.  (32)  are  of 
a  class  for  which  asymptotic  approximations  may  be  obtained 
for  either  large  (t]  ->  w)  or  small  (t)-^0,  f-^0)  values  of  rj;  As 
a  result  of  the  integrability  of  the  functions  g(cD)  the 
major  contribution  to  the  integrals  (32)  for  large  values  of 
T]  comes  from  the  vicinity  of  the  origin,  and  Laplace's  method 
as  demonstrated  in  Ref.  [12]  is  applicable. 


Case  a 

In  this  case  the  contributions  from  both  potentials  are 

of  the  general  form  (32),  I  =  I  +  I  .  Consider  the  typical 

1  2 

contribution  from  the  equivoluminal  potential  ^  .  In  this 
case 


-16 


and 


h,  (co)  = 

+  \  T) 


_  i  i  ^  -  i  i)  + 


M*^(l  -  m) 

T  -V'^« -  i.H.0(u.=  ).(34) 


Following  the  procedure  outlined  in  [12]  using  only  the 
first  term  of  the  expressions  (33)  ^rid  (3^)  o^ie  finds  for 


V  ^  c 


o  T]a3h_(cD)  N 

I  '  *11 

=  I  g_(cD)e  d(j3  ~  2^ 


n=o  (tip;l  -  iO 


n+  1 


(35) 


where 


Tia3[h_(co)-h_(0)  ] 


g_(cD)e 


=  E 

n=0 


A^co 


(36) 


Similarly^ 


«>  T)CDh  ( (Jo) 


g_l_(cD)e 


n!B 


N 

do)  ~  E  — 

n=0  (iiP*  +  i|) 


n 


n  +1 


(37) 


where 


g_j_(a))e 


•ncD  [h^((JD)-h_^(0)] 


E  B  CD^ 

n=0 


(38) 


Substitution  of  the  appropriate  expansion  coefficients 
and  into  Eqs.  (35)  (37)  leads  to  expressions  for 

the  stresses  and  accelerations.  Retaining  only  the  leading 
non-vanishing  coefficient,  A^  or  B^^  ,  respectively,  one 
obtains  for  V  c* 
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and 


Ti(B[h  (a3)-h  (0)]  B  ” 

g^(ai)e  ^  +  S  B^o) 

n=U 


respectively.  The  expressions  resulting  from  the  terms  other 
than  A  ^  and  B  ^  were  previously  obtained.  Only  the  in¬ 
tegrals 


u 

=/ 


dcD  + 


dco  (42) 


=  / 


o  ^  tiP^\  -ii 


—  le  -  e  Je  dw  -  I  —ie 


■Tip*co  -tiP*u5\  -i^cu 


need  be  considered. 

By  combining  the  two  integrals  in  I  one  finds 

3 


r 

=  2i 


-T)P*(I) 

00  I'^Jj 


e _  sin  ten 


dcD  =  2i  tan' 


In  a  similar  manner  one  finds 


I  = 

4 


-T)P*cD  -Tlpjejox 


-2  Me  -  e 


costen 


-log  .  (45) 


Substitution  of  the  appropriate  coefficients  A  or  B 

n  n 

into  Eqs.  (35) j  (37)j  (44)  and  (45)  leads  to  expressions  for 


the  stresses  and  accelerations.  Retaining  only  the  leading 


non-vanishing  coefficient,  or  ,  respectively,  one 

obtains  for  V  =  c* 

H. 


(2  -Mj2) 

^xx  ~  27TVTA* 

R. 


(2+M^2_  2M£2)tan-i  -  ( 2  -  )  tan"^  ^"1 


yy 


-  2M^^)(2  -  ^  ^ 

2'nVTAS  B*z 

K  L* 


(2  - 

^zz  27rV'IA* 

R 


P£(2- 


tan  ^  -  tan~^ 

P*z  pjz 


x^  + 


] 


a.._ - —  log 


xz  27A7TA* 

R 


x'^  +  P?‘='z' 
J_l 


(46) 


u  ~ 


27rmtj,'IA* 

R 


/  O  _  M*  ^  ^  ^  „  O  A*  2  ^ 

j  3C^  +  P*^z2  x^  + 


w 


27Tmjj,TA* 

R 


Xi 


p: 


*272 

T  ^  J 


where 


m: 


—  ^ 
n  ’  ^  m  ’  Pl 


1  -  =  1  -  , 


and 


*  _  M*a  1.  ~.  ^ 


A*  = 


m 


(2  - 


..*p  1  -  n  ,  ..*p  1  -  m 

"l  i^+»T  Spp-_| 
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The  result  (46)  differs  from  Eqs .  (39)  in  two  major  re¬ 
spects.  In  the  latter,  stress  and  acceleration  decay  as 
and  ^  respectively,  where  r^  =  ,  while  if 

V  =  c*  ,  the  stresses  do  not  decay  at  all  with  respect  to  r 
and  the  accelerations  decay  as  .  The  situation  is  some¬ 

what  comparable  to  the  resonance  of  a  damped  oscillator. 

The  second  significant  difference  between  Eqs.  (46)  and 
(39)  concerrs  a  matter  of  symmetry.  Quantities  and 

u^  which  were,  respectively,  symmetric  and  antisymmetric 
with  respect  to  x  when  V  ^  c*  ,  Eq.  (39)j  are  antisymmet- 

Xv 

ric  and  symmetric,  respectively,  when  V  =  c*  . 

R 

It  is  noted  that  there  is  no  sudden  change  in  behavior 
concerning  decay  and  symmetry  depending  on  whether  V  =  c* 

R 

or  V  ^  c*  ,  but  there  is  a  gradual  transition.  This  transi- 
K 

tion  in  the  vicinity  of  V  =  c*  is  not  contained  in  the 
present  results.  It  would  require  retention  of  more  than 
the  first  significant  term  in  the  asymptotic  expansion. 

Case  b 

In  this  case  the  exponent  of  the  irrotational  potential 
cE)  is  identical  with  that  for  Case  a.  Therefore,  Eqs.  (35) 
and  (37)  are  used  again  to  obtain  the  irrotational  contribu¬ 
tions  to  the  stresses  and  accelerations.  The  exponent  of 
the  equivoluminal  potential  ¥  is,  however, 

h_(u3)  =  h_|_(co)  =  h(cD) 


=  -  i  = 

^T  T) 


■ 


=  i(i4 


l'l^^(l-m) 

2nim* 


+  0(fx>)‘ 
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Since  the  first  term  of  this  expression  is  purely  imagi- 
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In  a  similar 

manner  one  finds 

B  exp 

—  — - - - 

'.^K  -  i 


provided  0  ,  while 


B  <  - -  + 


T  mmj 


i(nm*  -  i) 

>  nM*^  — -~- 
■  *  T  2iiun^ 


( -  o 

'  T  imn* 
T. 


Erfc  - 


-  i) 


if  B^  =  0  ,  B^  0  . 

Substitution  of  the  appropriate  expansion  coefficients 
and  B^  into  Eqs.  (35),  (37)^  (50),  (51)^  (52)  and  (53) 
leads  to  expressions  for  the  stresses  and  accelerations. 

The  leading  terms  of  these  expressions  (except  for  a  )  are 
of  the  type 


(x-m*z)2 


(x-m^z)' 

c(x-m*z)  "'5  z 


Erfc 


i  V^(x-ra*z) 


where  q  is  a  positive  constant.  All  other  terms  contain, 
with  or  without  the  exponential  factor,  higher  negative 
powers  of  z  such  that  the  response  is  only  of  consequence 


where 


=  —  6*^  =  1  -  m*2  ^2  _  w*2  _  n 

A*  =  (2  -  M*^)^  ,  AJ  =  (2  -  M^2)2  +  Umjp*  , 


and  the  contributions  of  the  irrotational  potential  $  have 
been  neglected  in  all  expressions  except  those  for  . 

Case  c 


In  this  case,  since  h_(a3)  =  h_|_(cD)  and  g  (cd)  =  g_|_(a)), 
Eq.  (32)  may  be  written  in  the  form 


I 


■ntoh(u)) 

e 


dco  . 


00 


(55) 


Since  the  exponents  are  of  the  form  of  Eq.  (47) ^  the  second 
terms  of  their  expansions  must  be  retained.  Thus 


I 


exp 


ioD(TimJ  -  i) 


1  -  m 
2mm*  _ 


do) 


(56) 


where  the  are  defined  in  Eq.  (36) . 

If  the  order  of  summation  and  integration  are  inter¬ 
changed,  each  term  is  recognized  as  a  Fourier  transform 
whose  inverse  may  be  taken  from  the  tables  [13].  Therefore, 
one  finds  for  the  nth  term 


exp 


(nm^  -  i)' 


*2 


1  -  m 


mm, 


T  J 


(57) 
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Substitution  of  the  appropriate  expansion  coefficients 
into  Eq.  (57)  leads  to  expressions  for  the  stresses  and 
accelerations.  The  leading  term  in  these  expressions  is  of 
the  type 


where  q  is  a  positive  constant.  All  other  terms  contain 
higher  negative  powers  of  z  and  an  exponential  factor  of 
the  above  type  such  that  the  response  is  only  of  consequence 
for  ratios  ^  near  m*  or  m^  .  To  obtain  the  significant 
part  of  the  response  it  is  therefore  sufficient  to  retain 
only  the  first  non-vanishing  coefficient  A^  .  At  other  lo¬ 
cations  the  response  may  be  deemed  to  vanish  in  first  approx¬ 
imation.  The  expressions  for  the  stresses  and  accelerations 
thus  found  for  ^  Z  ir^  are 


(2  + 

2A*  ^/zttVTM^^ 


exp  - 


2nm£ 


(x  -  m^z)^ 


-  2M^2) 

(2  +  -  2M^^)  *^xx 


(58) 
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V. 


Asymptotic  Approximations  for  Small  Values  of  t]  and 


For  small  values  of  and  $  the  asymptotic  evaluation 

gives  in  first  approximation  the  rather  obvious  result  that 
the  stresses  are  equal  to  those  in  an  elastic  medium  having 
the  properties  of  the  unrelaxed  viscoelastic  material.  There- 
fore,  no  further  analysis  will  be  required  except  for  the  spe¬ 
cial  situation  V  =  c  ^  the  Rayleigh  velocity  defined  by 
A  =  (2-mJ)2  -  =  0. 

In  the  special  situation  V  =  the  functions  h(cD)  and 
g((X))  exp[T)cu(h((Xi)  -  h(0))]  are  expanded  in  powers  of  ^  ,  with 

CO 

the  result; 

h  (cu)  =  -  i  -^  +  O(-) 

h_^(aj)  =  -  p  -i  ^  +  O(-) 

■ncu[h_(a3)-h_(0)  ]  ^  (60) 

g(a))e  =AcD^  +  Aa)  +  A  +  2Ato^ 

2  ^  °  n=l 

Tiuj[h_j_(a))-h_j_(0)  ] 

g.(a))e  =BC3^+Boo  +  B+  ZB  O)""  . 

210  -n 


The  value  of  the  integrals,  in  this  case,  is  principally 
due  to  the  contribution  from  large  values  of  t),  such  that 

-N  T]aih_(a))  “  'rjaih_^((D) 

I  =  I  +  I  :::  J  g_(a))e  do)  +  J  g_^{a))e  dco  (6l) 

-00  N 


where  N  is  a  large  but  finite  constant.  Using  the  trans- 
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formation  =  s  and  retaining  only  the  leading  non-vanish 
ing  coefficient  or  of  the  appropriate  power  series 

expansion  one  obtains 

(62 

provided  7^  0,  while 

(63 

if  A  =  0,  A  ^0. 

2  1 

For  small  values  of  t]  the  limit  -t^N  may  be  replaced 
by  zero.  This  leads  to  the  result 


-i  / 


(Pt  - 


A  s  e 
1 


2a. 


(64, 


if  A  7^  0,  or 
2 


(l^T  -  i?) 


(65; 


if  A  =  0,  A  ^  0. 


The  error  introduced  by  changing  the  upper  limit  of  in¬ 
tegration  is  NM  where  M  is  the  maximum  value  of  the  func¬ 
tion  gjcn)  exp(Tia.[h_(co)-h_(0)])  in  the.  interval  0  >  cn  ^  -  k 


This  error  can  be  made  as  small  as  desired  relative  to  the 
value  of  I  by  selecting  sufficiently  small  values  of  ^ 

7 

and  Ti . 

The  integral  I  may  be  treated  in  a  similar  manner  with 
8 

the  result 

2b 

I - ^ -  (66) 


provided  B  ^  Oj  or 
2 

I 

a 


(67) 


if  B  =  0,  B  ^  0. 

2  1 

Substitution  of  the  appropriate  expansion  coefficients 
and  B^  into  Eqs .  (64),  (65) j  (66)  and  (67)  leads  to 
expressions  for  the  stresses  and  accelerations.  Retaining 


only  the  leading  non-vanishing  coefficient  one  obtains 
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u  ~ 


Jj 


TTiaA 


R 


z(p^z^-  3x^)  z(pjz^-  3x2) 

(2-  Mj)  — i - rr—  -  2P' - - - 


(x=+ 


(x2+  pjz2)3 


W  ~- 


TTM- 


x(3P£z2-  x2)  x(3PJz2-  x2) 


where 


and 


P?  =  1-M^ 


-  2 


2  _ 


r. 


^R  =  Mj(l-m)(2-Mj)-  Pj.p^ 


^^(l-n)  Mj(l-m) 


L 


T  -> 


(68) 

cont , 


The  situation  here  is  quite  similar  to  the  one  at  large 
distances  when  V  =  c*  discussed  in  connection  with  Eqs.  (46). 
At  small  distances  r  0,  if  v  7^  c^^  the  stress  and  accelera¬ 
tion  increase  as  ^  and  ,  respectively^  [Ref.  4],  where 
r^  =  while  if  v  =  c^  the  stress  and  acceleration 

increase  as  and  ,  respectively.  This  is  again  some¬ 

what  comparable  to  the  resonance  of  a  damped  oscillator. 

The  various  quantities  defined  by  Eqs.  (68)  show  symmetry 
or  antisymmetry  with  respect  to  x  opposite  to  the  behavior 
of  the  equivalent  quantities  in  the  elastic  solution  for  ^  ^ 
[Ref.  4]. 

As  in  the  case  V  =  c*,  there  is  no  sudden  change  in  be- 
havior  concerning  decay  and  symmetry^  depending  on  whether 
V  =  c„  or  V  7^  c  ,  but  there  is  a  gradual  transition.  This 
transition  in  the  vicinity  of  V  =  c^^  is  not  contained  in 
the  result  (68) .  It  would  require  the  retention  of  more  than 
the  first  significant  term  in  the  asymptotic  expansion. 
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VI.  Discussion  of  Results 


The  asymptotic  approximations  obtained  permit  a  discus¬ 
sion  of  the  differences  between  the  response  for  elastic  and 
viscoelastic  materials.  Deferring  discussion  of  the  special 
situation  V  =  c  ,  the  asymptotic  approximations  for  points 
at  small  distances  are  found  to  be  simply  the  elastic  sub¬ 
sonic  solutions,  "small"  or  "large"  to  be  interpreted  in  com¬ 
parison  with  the  products  of  the  relaxation  time  T  and  the 
typical  wave  velocities  c^  •>  ^  “^p  "^s  ‘  Using  the 

smallest  and  largest,  respectively,  the  asymptotic  approxi¬ 
mations  for  "small"  distances  apply  if 

r  «  c*T  ,  (69) 

while  the  asymptotic  approximations  for  "large"  distances 
apply  if 

r  »  CpT  .  (70) 

The  differences  between  the  elastic  and  viscoelastic 
situations  appear  in  the  far  field.  These  differences  de¬ 
pend  on  the  relative  values  of  the  velocity  V  and  the  wave 

velocities  c*  and  c*  .  Three  cases,  previously  designated 
p  s 

a  ,  b  and  c  ,  have  been  considered. 

Case  a  V  <(  c* 

In  this  case  the  velocity  of  the  load  is  less  than  the 
propagation  velocity  of  any  plane  wave  in  the  material,  and 
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the  effect  of  viscosity  on  the  stresses  and  accelerations  is 

not  radical.  Consider,  for  example,  the  stresses  as  func- 

tions  of  the  distance  r  =  ,  but  such  that  —  is 

constant.  Fig.  6.  In  the  elastic  material  the  stresses  are 

equal  to  —  ,  where  c  depends  only  on  the  value  of  —  . 

3r  z 

In  the  viscoelastic  material,  excluding  the  special  case 

C  “ 

V  =  c*  ,  the  stresses  may  still  be  written  ,  but  c  has 

a  slight  dependence  on  r  .  For  small  values  of  r  ,  c  =  c, 

while  for  large  values  of  r  ,  c  ,  Eqs,  (39) j  approaches  a 

slightly  different  value.  As  an  example,  the  stress 

along  —  =0  is  shown  in  Fig.  7  a  function  of  z  for 

the  specific  set  of  parameters  =0.3^  m=0.5j  ^“3*  • 

At  large  distances  the  ratio  ,  in  this  example,  approaches 

the  value  I.56.  The  value  of  this  ratio  depends  on  the  value 

of  —  and  may  be  negative.  Due  to  the  fact  that  only  the 
z 

asymptotic  values  have  been  obtained,  the  details  of  the  de¬ 
pendence  of  c  on  r  are  not  known . 

The  special  situation,  V  =  c*  ,  may  be  compared  to  the 

X\ 

response  of  a  damped  oscillator  at  resonance,  where  stresses 
and  displacements  are  much  larger  than  at  non-resonant  fre¬ 
quencies.  Similarly,  for  this  valuej  of  V  the  stresses  in 
the  far  field  are  much  larger  than  for  other  values  of  V  . 
They  no  longer  decay  as  ^  ,  but  approach  constant  values 
depending  only  on  the  ratio  ^  .  The  accelerations  which 
previously  decayed  as  now  decay  as  .  In  addition, 

quantities  previously  symmetric  with  respect  to  the  variable 
X  become  antis;innmetric  and  vice  versa,  as  discussed  in 
connection  with  Eq.  (46) . 
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Case  b 


c*  >  V  >  c* 

In  this  case,  where  the  velocity  V  exceeds  the  velocity 

of  shear  waves  for  low  frequencies,  ->  0  ,  the  far  field 

differs  radically  from  the  near  field,  and  from  the  far  field 

in  the  elastic  material.  Before  proceeding  further,  the 

reader  is  reminded  of  the  character  of  the  supersonic  elastic 

solution,  i.e.  when  V  )►  c  ]>  c  .  In  this  case  the  stresses 

p  s 

are  discontinuous .  They  vanish  everywhere  except  at  points 

where  —  has  one  of  the  values  defining  the  P  and  S 

fronts.  At  these  points  the  stresses  become  infinite^.  In 

the  elastic  trans-sonic  case,  the  stresses  are  the  sum  of  a 

smooth  function  —  c (— j  and  of  one  which  is  discontinuous 

r  \z/ 

at  the  S  front. 

In  the  viscoelastic  case  under  consideration  the  response 
in  the  near  field  is  again  given  by  the  subsonic  elastic  so¬ 
lution,  just  as  in  Case  a.  The  stresses  are  of  the  form  — 
where  c  is  a  smooth  function  of  —  .  In  the  far  field,  how- 
ever,  the  stresses  vanish  everywhere,  in  first  approximation, 
except  in  the  vicinity  of  —  =  m*  .  The  presence  of  the 

Z  J. 

viscous  effects  changes  the  response  in  the  far  field  from  a 
smooth  one,  into  one  somewhat  similar  to  the  one  in  the  elas¬ 
tic  trans-sonic  case.  In  the  viscoelastic  material  large  yet 
finite  stresses  occur  within  a  small  angle  in  the  vicinity 
of  a  critical  direction  defined  by  —  =  m*  .  In  essence,  the 

^  While  the  stress  becomes  infinite,  the  integral  of  the 
stress  across  either  front  is  finite  and  independent  of  r 
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viscoelastic  material  focuses  the  response  in  this  critical 
direction.  The  focusing  is  not  sharp^  but  decays  exponentially 
as  shown  in  Eqs.  (5^).  Along  the  critical  direction  the 


stresses  j  except  o  ,  are  of  the  form  -7= 

yy  V  z 


As  an  example^ 


the  stress  a  along  —  =  m*  is  shown  in  Fig.  8  as  a  func- 

ZZ  Z  T 

tion  of  z  for  the  specific  set  of  parameters  =  0.75  :> 

m  =  0.5  ,  K  =  ^  . 


Case  c 


V  >  c*  >  c* 
P  ^  s 


In  this  case  the  far  field  response  in  the  viscoelastic 
material  again  differs  radically  from  that  in  an  elastic  ma¬ 
terial.  The  situation  is  somewhat  similar  to  that  in  Case  b^ 

but  there  is  focusing  in  the  two  directions  —  =  m*  and 

z  Lj 

X 

—  =  .  It  is  to  be  noted  that  Case  c  is  only  possible  if 

K  )>  27 ^  ;  this  can  occur  only  for  a  very  small  or 

negative  value  of  Poisson's  ratio  in  combination  with  m  <(<(l. 


Special  Case  V  =  c  ,  the  Velocity  of  Rayleigh  Waves 

_ _  X\ 


In  this  case  the  elastic  solution  does  not  exist  at  all. 

The  stresses  increase  in  magnitude  as  V  approaches  c  , 

just  as  in  the  resonance  of  an  undamped  oscillator.  The 

presence  of  viscous  effects  leads  to  a  finite  response,  Eqs. 

(68).  The  behavior  of  the  solution  for  small  values  of  r 

differs  radically  from  the  elastic  solution  for  V  7^  c  in 

R 

two  respects.  Stresses  and  accelerations  in  the  latter  case 
increase  for  r  0  as  ^  and  ,  respectively,  while 

Eqs.  (68)  indicate  an  increase  as  and  ,  respectively. 
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In  the  vicinity  of  r  =  0  the  response  for  ^ 
therefore  an  order  of  magnitude  larger  than  for  V 
In  addition j  there  is  a  change  from  symmetric  to  antisymmet¬ 
ric  behavior  with  respect  to  the  variable  x  and  vice  versa. 
The  changes  are  quite  similar  to  those  encountered  in  the 
far  field  when  V  =  c*  , 

The  far  field  for  V  =  c  can  be  obtained,  depending 

on  the  material  properties,  from  Case  a,  b,  or  c.  The  fact 

that  V  =  c  does  not  create  a  special  situation  in  the  far 
R 

field,  r  ^  00  . 
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FIG.  3  BRANCH  POINTS  AND  CUTS  USED  IN  CASE  a 


FIG,  4  BRANCH  POINTS  AND  CUTS  USED  IN  CASE  b 
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